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ABSTRACT 

Context. Cosmic rays are present in almost all phases of the ISM. PAHs and cosmic rays represent an abundant and ubiquitous 
component of the interstellar medium. However, the interaction between them has never before been fully investigated. 
Aims. To study the effects of cosmic ray ion (H, He, CNO and Fe-Co-Ni) and electron bombardment of PAHs in galactic and 
extragalactic environments. 

Methods. We calculate the nuclear and electronic interactions for collisions between PAHs and cosmic ray ions and electrons with 
energies between 5 MeV/nucleon and 10 GeV, above the threshold for carbon atom loss, in normal galaxies, starburst galaxies and 
cooling flow galaxy clusters. 

Results. The timescale for PAH destruction by cosmic ray ions depends on the electronic excitation energy Eq and on the amount of 
energy available for dissociation. Small PAHs are destroyed faster, with He and the CNO group being the more effective projectiles. 
For electron collisions, the lifetime is independent of the PAH size and varies with the threshold energy Tq. 

Conclusions. Cosmic rays process the PAHs in diffuse clouds, where the destruction due to interstellar shocks is less efficient. In the 
hot gas filling galactic halos, outflows of starburst galaxies and intra-cluster medium, PAH destruction is dominated by collisions with 
thermal ions and electrons, but this mechanism is ineffective if the molecules are in denser cloudlets and isolated from the hot gas. 
Cosmic rays can access the denser clouds and together with X-rays will set the lifetime of those protected PAHs. This limits the use 
of PAHs as a 'dye' for tracing the presence of cold entrained material. 

Key words, dust, extinction - cosmic rays - Galaxies: halos - Galaxies: starburst - ISM: jets and outflows - shock waves - Galaxies: 
cooling flows - intergalactic medium 
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0^ 1. Introduction Unfortunately the lack of theoretical studies on PAH pro- 

cessins in shocked regions combined with the difficulty in dis- 

^ A rernarkable characteristic of the IR emission features is that ^^e PAH features intrinsic to the shocked region with 

^ they dominate the mid-lR spectrum of alinost all objects associ- ^j^^^^ ^^^^ surrounding material makes the interpre- 

^ ated with dust and gas and illuminated by UV photons, rang- ^^^-^^ of such observations ra ther complicated. In our previous 

O ; . ^.^^''T.^o'! '° ultraluminous infrared galaxies (see ^^^j^^ jMicelotta et alJl20lo3bl hereafter MJTa and MJTb) we 

iTieMsl I20Q81 for a recent review). At present, these feature ^^^^.^^ ^^^^j^^j ^^^^^.^ molecules in interstellar shocks 

are (alinost) universally attributed to the IR fluorescence of far- ^^^j^ velocities between 50 and 200 km s^^ and in a hot post- 

ultraviolet (FUV)-pumped Polycyc ic Aromatic HydrocaAon ^j^^^,^ ^^^j^ th e Herbig- Haro je ts in the Orion and Vela 

(^) molecules containing 50 - 100 carbon atoms dTieleiii ^^^^ fo rming regions (iPodio et alJl2006h . in the local interstellar 

PAHs have recently been detected m association with ^j^^^ (Slavin 20oi and in flie outflow of the starburst galaxy 

shocked hot gas, but it is difficul t toe|ablish a clear comec- (Engelbracht eTdlM) . We found that interstellar PAHs 

tion between the two. | Tappeetal. | (l2006|) detected spectral fea- ^ ^^^^.^^ ^^^^^^ ^.^^ velocities greater than 

tures in the emission of the supernova reinnant N132D m the joOkms' and larger PAHs (A^c = 200) are destroyed for shocks 

Large Mage lanic Cloud ,^h they attribute to emission by ^.^^ velocities > 125 km s '. Even where destruction is not 

large PAHs. [ Reach et al. | have identified four supernova ^^mplete, the PAH structure is likely to be severely denatured 

remnant with R£^rs maybe indicating PAH einission, and ^ ^j^^ j^^^ important fraction (20 - 40%) of the carbon 

l Aiidersenet al. 1 (l2007|) investigated the presence of PAHs in a ^^^^^ .^^1 p^j^ ^-^^^-^^^ ^^^^^ ^ ^ 

subset of galactic supernova remnants in the GLIMPSE survey. ^j^^ ^ ^^^^^^^ ^ q ^^_3^ ^ 

PAHs have also been observed interwoven with the X-ray emis- ^q7 k), typical of the coronal gas in galactic outflows, PAHs are 

sion a rising from the bi polar outflow of the starburst galaxy i^^j „ destroyed by electron collisions, with lifetimes mea- 

M82 (| Armusetal | |2007^ and in the high-latitude coronal gas sured in thousands of years, i.e. orders of magnitude shorter than 

?f the edge-on ^^I^.NGC 5907 NGC 5529 and NGC 891 ^^e typical lifetime of such objects. 
dlrwin & Maddea2006i:llrwin et al.ll2007HWhalev et al1l2009l) . 

Cosmic rays (CRs) are an important component of the ISM, 

Send offprint requests to: E.R.MicelottR contributing considerably to its energy density (^ 2 eV cm"-*. 
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lTielenj2005HPadovani et al.ll2009l) . CRs consist mainly of rela- 
tivistic protons, a-particles (~ 10%), and heavier ions and elec- 
trons (~ 1%). The spectrum (intensity as a function of the en- 
ergy) of the ionic CR component measured near the Earth spans 
from ~ 100 MeV to ~ 10~" eV, and decreases steeply with en- 
ergy. The spectrum of the electronic compo nent is even steeper 
and ranges from ~ 600 MeV to 10^ GeV (lip & Axfordlll985t 
[Gais ser & Stanev 2006) . 

The lowest-energy CRs in the ISM, with energy between 
5 MeV and few GeV, are excluded from the heliosphere or 
severely slowed down by the solar wind. Hence, the y cannot 
be dir ectly observed even with far-ranging spacecraft ( Shapirol 
1 19961) and th eir spectra have to be evaluated theoretically 
(IShapird[T99Th . 

CRs with energy up the few 10'^ eV (the "knee" observed 
in the spectrum: Drurv 1994) are thought to be produced in the 
Galaxy, mainly by supernova shocks in the disk. Because of their 
charge, CRs are tied to the galactic magnetic field and are con- 
fined to a spheroidal volu me with radius of ~ 20 kpc and half - 
thickness of ~ 1 - 15 kpc ( lGinzburg|ll9"8llShibata et alj|2007h . 
with a small but finite escape probability. The magnetic field ran- 
domizes the trajectories of CRs as they propagate through the 
Galaxy, so their distribution is almost isotropic except close to 
the sources. 

From the point of view of PAH destruction, CRs have then 
two interesting characteristics; first, for energies up to 10 GeV 
they can efficiently transfer energy to the PAH, with possible 
consequent destruction (see Sects. 2 and 3); second, they per- 
meate almost homogeneously the ISM and can penetrate into 
regions such as dense clouds which are otherwise not much af- 
fected by high temperature ions and electrons (MJTb). 

The aim of this work is to quantify the destructive potential 
of CRs and to compare it with other mechanisms (interstellar 
shocks, collisions within a hot gas. X-ray and FUV absorption), 
in galactic and extragalactic environment. 

The paper is organized as follows: Sect. 2 and Sect. 3 de- 
scribe the treatment of high energy ion and electron interactions 
with PAHs, Sect. 4 presents the CR spectra adopted for our study 
and Sect. 5 illustrates the calculation of the collision rate be- 
tween PAHs and CRs. We present our results on PAH destruction 
and lifetime in Sect. 6 and discuss the astrophysical implications 
in Sect. 7, summarizing our conclusions in Sect. 8. 



2. High energy ion interactions with solids 

2.1. Collisions with high energy ions 

To describe the effects of high-energy ion collisions with PAH 
molecules we adopt a similar approach to that used in our previ- 
ous work (MJTa and MJTb), based on the theory of ion interac- 
tion in solids. Ions colliding with a PAH will excite the molecule 
and the excitation can lead to fragmentation or, alternatively, the 
excess energy can be radiated away. A calculation of the frag- 
mentation process thus consists of two steps: 1) the calculation 
of the excitation energy after collision, 2) the probability of dis- 
sociation of an excited PAH. The former is discussed here. The 
latter is described in section Sect. 2.2. 

Energetic ion s can interact with a s ohd through two simulta- 
neous processes dLindhard et al.l[T963h : the nuclear stopping or 
elastic energy loss, where the energy is directly transferred from 
the projectile ion to a target nucleus via a binary elastic collision, 
and the electronic stopping or inelastic energy loss, consisting of 
the energy loss to the electron plasma of the medium. For high 
energies characteristic of CRs (above ~ 1 MeV/nucleon) we only 



need to be concerned by electronic stopping, whi ch can be well 
described by the Bethe-Bloch equation (see e.g. IZieglerl[l999l 
and references therein). 

The Bethe-Bloch equation was derived considering the elec- 
tromagnetic interaction of an energetic particle with the electron 
plasma of a solid. A PAH molecule is not a solid but its large 
number of delocalized electrons can be treated as an electron gas 
(see MJTb and references therein). It is therefore appropriate to 
consider the energy loss to such an electron plasma in terms of 
the Bethe-Bloch electronic stopping power, S (energy loss per 
unit length). When S is known, a specific procedure has to be 
applied to calculate the effective amount of energy transferred 
into a single molecule, taking into account the finite geometry of 
the PAH (see Sect. 3). 

The stopping of high velocity ions in matter has been a sub- 
ject of interest for more th an a century, starting with the work of 
Marie Curie in 1898-1899 ( ICurielll900h. For a the oretical review 
on the topic, we refer the reader to lZieg ie3 (ll999h and references 
therein. In the following, we summarize the basic methods and 
equations for evaluating the stopping power of high energy ions 
and we illustrate the modifications introduced into the theory in 
order to treat the interaction with PAH molecules. 

For high energy light ions (H, He and Li above 1 
MeV/nucleon), the fundamental relation describing the stopping 
power, S , in solids is the relativistic Bethe-Bloch equation, com- 
monly expressed as 



S = -^ZfL(J3) 



kZ2 2 



Z\ [Lo(y6) + ZyUifi) + Z\L2{Si)..\ 



where 



K = \nr\m^(?', ro = e^lmsC^ (Bohr electron radius) 



(1) 
(2) 

(3) 



ni^ is the electron mass, Z] and Z2 are the projectile and target 
atomic numbers respectively and - v/c is the relative projectile 
velocity. The term LiJJ) is defined as the stopping number, and 
its expansion in Eq. 2 contains all the corrections to the basic 
ion-electron energy loss process. 

Th e first term Lp includes the fundam ental Bethe-Bloch re- 
lation (lBethelll930l 119321: lBlochlll933allbl) for the stopping of 
high-energy ions, tog ether with th e two corrective terms C/Z2 
and 512 introduced bv lFanol(ll963h 



Lo^m-^ - ln</) - ^ 



where 



m = In 



1 -y82 



(4) 



(5) 



C/Z2 is the shell correction and takes into account the fact 
that as soon as the projectile loses energy into the target, its ve- 
locity decreases from relativistic values, thus the Bethe-Bloch 
theory requirement of having particles with velocity far greater 
than the velocity of the bound electron is no longer satisfied. 
In this case, a detailed accounting of the projectile's interaction 
with each electronic orbital is required. The shell correction is 
important for protons in the energy range of 1 - 100 MeV, with 
a maximum contribution of about 10%. 

ln(/), which is part of the original Bethe-Bloch relativistic 
stopping formula, represents the mean ionization, and corrects 
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for the fact that the energy levels of the target electrons are quan- 
tized and not continuous. 

(5/2 represents the density effect and provides the correction 
to the reduction of the stopping power due to polarization ef- 
fects in the target. The dielectric polarization of the target mate- 
rial reduces the particle electromagnetic field from its free-space 
values, resulting in a variation of the energy loss. The density 
effect becomes important only for particles with kinetic energies 
exceeding their rest mass (938 MeV for proton). 

The term Zi Li takes into account the Barkas effect. This is 
due to the target electrons, which respond to the approaching 
particles slightly changing their orbits before the energy loss in- 
teraction can occur. 

The term ZjLa is the Bloch correction and provides the 
transition between the two approaches used to evaluate the en- 
ergy loss of high-energy particles to ta rget electrons: the clas- 
sical Bohr impact-parameter approach (iBohjl 19131 Il915h . and 
the quantum-mech anical Bethe momentum transfer approach 
(lBethelll930lll932h . 

Both the Barkas and Bloch corrections are usually quite 
small and contribute less than few percent to the stopping at en- 
ergies from 1 MeV/nucleon to 10^ MeV. 

The term Lq can be evaluated as a function of the energy of 
the incoming ion using the following definition (IZiegledll999h : 



1 



1 



(6) 



yc/ [1 -H £(GeV)/ 0.931494 Mi(amu)]' 

Using Eq. |4] we can rewrite the expression for the stopping 
power (Eq. 2) in the following way 



m - ^ 

Z2 



ln</) - - 



Zi Li ( Barkas ) -i- ZjL2{ Bloch )} (7) 



9 



For the stopping power in units of eV/(10'^ atoms/cm^) 
the pref actor constant has the value k = 5.099x10 while for 
the stopping units of eV/A the above prefactor has to be mul- 
tiplied by N/IQ^^, where is the atomic density of the target 
(atoms/cm^). In our calculation we adopt the values for amor- 
phous carbon (Z2 = 6, M2 = 12, A^ = 1.1296x10^^^ atoms/cm^^), 
which yields k = S.TSOSxlO"'^. 

An empi rical expressi on for the Barkas correction term 
is available (IZiegleii(l999l) . whereas the Bloch correction can 
be evaluated usi ng the Bichsel parametrisation (lBichsellll990h . 
lBonderupl(ll967h estimates that (/) = 1 1 .4 Z2 (eV), but this is not 
always in agreement with experimental data and unfortunately 
there are no simple algorithms for the shell correction and den- 
sity effect. As a consequence, no simple analytical expressions 
for the stopping power are available. 

The SRIM program (Ziegl eret al.lll985h calculates accurate 
stopping powers from Eq. |7] using different methods to evalu- 
ate the corrective terms. The shell correction C/Z2 is the average 
of the values obtained from the Local Density Approximation 
theory (LDA) and Hydrogenic Wave Function (HWF) approach 
(Ziegler 1995^. The first is an ab initio calculation based on 
realistic solid state charge distributions, while the second uses 
parameterized functions based on experimental stoppin g data. 
For th e density correction 6/2, the values tabulated in lICRUl 
(1 19841) are used, while the term ln(/) is derived by adjust- 
ing the theoretical value obtained from the LDA theory (e.g. 
iLindhard & ScharfalT952l) in order to fit the sum ln(/) -1- C/Z2 
evaluated empirically from experimental stopping data. 



As already mentioned, Eq. |7] is valid for light ions, H, He 
and Li. The stopping power for ions with Zi > 3 is usually cal- 
cu lated using t he hea vy-ion scaling rule, as reported for example 
bv lKatzetal.l(ll972h 



s{Zi,/3) = s(p,p)[zi/z;f 



(8) 



where S{p, /3) is the stopping power of a proton at the same 
speed as the ion of atomic number Zi (Eq. |7]i, Z* and Z* are the 
effective charge numbers of ion Z[ and of a pro ton respectively, 
with the expression given by Barkas (see e.g. iHenriksen et al.l 
1I970I) . 



Z* =Zi [1 - exp(-125ySZj"^^^^)] 



(9) 



To verify the effective importance of the corrective terms to 
the Bethe-Bloch equation in our specific case of interest (H, He, 
C and Fe impacting on carbon), we compared the output from 
SRIM with the stopping power calculated from the following 
approximate equation 



S-'-^Zl [m - ln</)] 



(10) 



where we adopt the carbon mean ionization energy calculated by 
SRUVQ using the method described above, (/) = 79.1 eV. For the 
heavier ions C and Fe we use the proton stopping power from 
Eq.[T0]intoEq.[8] 

We find a maximum discrepancy between the two curves of 
~10%, very small with respect to the other astrophysical uncer- 
tainties and indicating that our approximate equation is adequate 
to describe the energy loss. In this specific case, not only the 
Barkas and Bloch corrections, but also the shell correction and 
the density effect play a marginal role. Fig.[T]shows the compari- 
son between the 'accurate' SRIM curve and the approximate one 
for H and Fe, our lightest and heaviest projectiles respectively. 

Using Eq.|8]and[T0]we can then calculate the electronic stop- 
ping power dE/dx of energetic ions and from this the energy loss 
to the PAH molecule (Sect. 3). Since the stopping power is a de- 
creasing function of the ion energy and increases quadratically 
with Zi, the major contribution to the energy loss will come from 
the less energetic particles for a given ion, and from the heavier 
species for a given velocity (Fig.[TJ. 

2.2. Ion energy loss and dissociation probability 

To calculate the energy transferred to a PAH during collisions 
with high energy ions we adopt the configuration shown in Fig.|2] 
(see also MJTb). 

The molecule is modeled as a thick disk with radius R given 
by the usual expression for the radius of a PAH: apah - R = 
0.9 V(^c) A^_wherejyc is the number of carbon atoms in the 
molecule (iOmont'1986'). For a 50 C-atom PAH, R = 6.36 A. The 
thickness of the disk, d ~ 4.31 A is the thickness of the electron 
density assumed for the PAH (see MJTb). 

The path I, through the PAH, along which the incoming ion 
loses its energy is defined by the impact angle § and by the di- 
mensions of the molecule. Inspection of Fig. |2] shows that, if 
I tan()9-)| < tan(a), = d/\ cosi?|, otherwise /(i?) = 2R/\ sm&\. 
The rigorous method to calculate the energy loss along l{§) 
takes into account the progressive slowing down of the projec- 
tile in traversing the target. This is explained in Sect. 3 of MJTb. 
Nevertheless, in the present case we consider high energy parti- 
cles for which the energy loss along the path I is small compared 
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Fig. 1. Stopping power of hydrogen and iron impacting on car- 
bon, as a function of the energy per nucleon of the ion. The lower 
validity limit of the Bethe-Bloch equation is 1 MeV/nucleon. We 
compare the output from the SRIM code, which computes all 
corrections, and our approximate equation which includes only 
the mean ionization correction. 




Fig. 2. Coordinate system adopted to calculate the energy trans- 
ferred to a PAH via electronic excitation by ion collisions. The 
molecule is modeled as a disk with radius R and thickness d. The 
trajectory of the incoming particle is identified by the angle &, 
while the angle a corresponds to the diagonal of the disk. 



to the initial energy (cf. Fig. [TJ. In other words, the energy of 
the incoming ion remains almost constant during the interaction, 
thus the amount of energy lost after travelUng the distance /()?) 
can be simplified to 



dE 

T{§) = — X 
dx 



(11) 



where dE/dx = S{&) is the stopping power (energy lost per unit 
length) from Eq. [10] and/or Eq. [8] and dE is the amount of en- 
ergy lost after travelling the distance d.Y. The difference between 
the results from the two methods is of the order of few percent, 
which fully justifies the use of approximation in Eq. [TT] 

The ion collision will leave the molecule electronically ex- 
cited. Internal conversion and/or intersystem crossing will trans- 
fer this excitation (largely) to the vibrational manifold. De- 
excitation can occur through two competing decay channels: 
emission of infrared photons and dissociation and loss of a C2 
fragment. The latter is the process that we are interested in be- 
cause it leads to PAH fragmentation. The emission of a C2 frag- 
ment is suggested by experiments on fullerene Ceo which have 
shown that the ejection of C2 groups is one of the preferred frag- 
mentation channel. Moreover, the loss of acetylene groups C2H2 
has been observed in small PAHs. In fact, in a PAH molecule, a 
side group C2H,, (with n - 0, 1,2) is easier to remove because 
only two single bonds have to be broken, while the ejection of a 
single external C-atom requires one single and one double bond 
to be broken, and for an inner C-atom from the skeleton three 
bonds needs to be broken. To quantify the PAH destruction due 
to ion collisions we need to determine the probability of dissoci- 
ation, P, rather than IR emission. For the detailed calculation and 
a discussion of the dissociation probability, we refer the reader to 
Sect. 4. 1 in MJTb. For the sake of clarity, the basic equations are 
reproduced here. The total dissociation probability is calculated 
by combining the rates for fragmentation (A:o exp[-£'o/^ 7"av]) 
and IR decay (A:iR/(«max + 1)) into the following expression 



^0 exp [-Eo/k Ta. 



[^IR/(nmax + 1)] + ^0 s\p [-E^ / k T^y] 



(12) 



where ^0 and £0 are the Arrhenius pre-exponential factor and en- 
ergy describing the fragmentation process respectively, A:ir and 
«max are the IR photon emission rate and number of IR photons 
(MJTb), and k is the Boltzmann constant. The temperature Tav is 
chosen as the geometrical mean between two specific effective 
temperatures of the PAH 



(13) 



In the microcanonical description of a PAH, the temperature, T, 
describing the excitation (for fragmentation purposes) is related 
to the internal energy, E', by 



T ^ 2000 



£'(eV) 



1 - 0.2 



£Q(eV) 
£'(eV) 



(14) 



where £0 is the binding energy of the fragment (lTielensll2005K 
The temperatures To and T„^^^ in Eq. [13] are the temperatures 
when the internal energy equals the initial transferred energy 
(£' = E) and when the internal energy equals the energy af- 
ter emission of «niax photons {E' - E- nmax x As) with Ae being 
the average energy of the emitted IR photon. For the number of 
photons, Hmax, required to be emitted to have the probability per 
step drop by an order of magnitude, we adopt 10, the average 
photon energy is set equal to 0.16 eV, corresponding to a typical 
CC mode, and the pre-exponential is set equal to 1.4 x 10'^ s ' 
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(MJTb). For the photo n emission rate we ad opt the typical value 
/tiR = 100 photons s ' dJochims et al.lll994 . 

From the above equations one can see that P depends on the 
binding energy of the fragment, £0, on the PAH size, A^c, and 
on the energy transferred, T, which in turns depends on the ini- 
tial energy of the projectile. For a fixed value of the transferred 
energy, the dissociation probability decreases for increasing £0 
and A^c because either more energy is required in the bond that 
has to be broken or because the energy is spread over more vi- 
brational modes and hence the internal excitation temperature is 
lower. On the other hand, the more energy that is deposited in 
the PAH, the higher is the dissociation probability. 

The fragment binding energy £0, which is a crucial param- 
eter in the evaluation of the dissociation probability is, unfortu- 
nately, presently not well constrained. As in our previous work 
(MJTb) we investigated the impact on the PAH destruction pro- 
cess for £0 = 3.65, 5.6 and 4.58 eV, the latter value is consistent 
with extrapolations to interstellar conditions and is our reference 
value. 

As mentioned at the end of Sect. 2. 1, the stopping power in- 
creases with the atomic number of the projectile and decreases 
with its energy. Thus, for a given pathlength, the energy trans- 
ferred, and therefore the dissociation probability, will be higher 
for low energy heavy particles. 




LABORATORY COORDINATES 

Fig. 3. Scattering geometry for the elastic collision of an elec- 
tron (mass me, impact parameter p, initial velocity v) on a mas- 
sive particle (mass M2, initial velocity zero). After the impact, 
the target particle is knocked out and recoils at an angle with 
respect of the initial direction of motion of the projectile elec- 
tron. 



3. Collisions with high energy electrons 

To model the interaction of high-energy electrons with PAH 
molecules, we refer to the formalism used to describe the irra- 
diatio n effects in sol id materials, in particular carbon nanostruc- 
tures (lBanhartlll999l) . In the collisions of high-energy (relativis- 
tic) electrons with nuclei, the screening effect of the surrounding 
electrons is negligible. The electron-nucleus interaction can thus 
be treated in terms of a binary collision using a simple Coulomb 
potential, applying the appropriate re lativistic corrections (e.g. 
iReimer & Braunlll989l:lBanhartlll999l) . 

If the energy transferred to the nucleus exceeds the displace- 
ment energy T^, i.e. the minimum energy required to produce a 
vacancy- interstitial pair which does not spontaneously recom- 
bine, the atom will be knocked out. If its energy is above the 
threshold value for further displacements, it can remove other 
atoms in its environment generating a collision cascade. 

In this description the "bulk" nature of the target enters only 
after the first interaction, when projectile and displaced atom 
propagate into the solid. Therefore, if we limit ourselves to the 
first interaction only, this approach can be applied to electron- 
PAH collisions and allows us to take into account the "molecu- 
lar" nature of our target. In fact, this is the same binary collision 
approach used to describe the nuclear interaction (elastic energy 
loss) in collisions between PAHs and relatively low energy ions 
in interstellar shocks (MJTa). 

A PAH is a planar molecule with tens to hundreds of car- 
bon atom. In this case the target nucleus is a single carbon in 
the PAH. If the energy transferred exceeds a threshold value the 
target nucleus will be ejected from the molecule. The displace- 
ment energy then has to be replaced by the threshold energy Tq, 
which represents the minimum energy to be transferred in or- 
der to knock-out a carbon atom from a PAH. The interaction 
between a high energy CR electron and a PAH occurs between 
the impinging electron and one single target carbon atom in the 
molecule. The energy is transferred from the projectile electron 
to a target carbon atom through a binary collision. Thus, each 
electron collision implies the loss of one single carbon atom, in 



contrast to ion collisions where each interaction causes the ejec- 
tion of a C2 fragment from the PAH molecule. The scattering 
geometry is shown in Fig. [3] After the collision, the target nu- 
cleus is knocked out and recoils at an angle with respect of the 
initial direction of motion of the projectile electron. The energy 
T transferred to the nucleus depends on the scattering angle 

r(0) = r„,axcos2© (15) 

The term r„,ax is the maximum transferable energy, correspond- 
ing to a h ead-on collision (0 = 0), and is given by the following 
equation dSimmon il965h 

2£ (£ + 2OTec2) 



where £ is the electron kinetic energy, nig is the electron mass 
and M2 is the target atomic mass. 

The total displacement/removal cross section, cr, i.e. the 
cross section for collisions able to transfer more than the thresh- 
old energy To, is defined as the integral over the solid angle of 
the differential cross section dcr/dO, which provides the proba- 
bility for atomic recoil into the solid angle dQ 

r®™- da 

cr= — 27rsin0d0 (17) 

Jo dn 

where = corresponds to the transfer of T,^^^ and = ©max is 
the recoil angle corresponding to the transfer of the minimun en- 
ergy Tq. The calculation of the total cross section cr for atom dis- 
placement/removal would require the analytical treat ment of the 
Mott scattering of a relativistic electron by a nucleus dMotll 19291 
Il932h . T he corresponding equations ha ve to be solved numeri- 
cally, but lMcKinlev & FeshbachI (Il948h found an analytical ap- 
proximation which provides reliable values of cr for light target 
elements such as carbon (under the assumption of an isotropic 
displacement/threshold energy). We adopt the formulation of the 
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analytical expression reported bv lBanhari (Il999h 



To 



'- max 1 2 



1 + 2na/3 



T 



1/2 



1 + 27Ta/3 + {p~ + nap) In 



(18) 



where Z2 is the atomic number of the displaced atom (in our 
case, carbon), is the Rydberg energy (13.6 eV), m^c^ is the 
electron rest mass (0.5 1 1 MeV), qq is the Bohr radius of the 
hydrogen atom (5.3x10"" m), yS = v/c, v being the velocity 
of the incident electron, and a - Z2/I37, where 1/137 is the fine 
structure constant. 

The term To is the minimum energy which has to be trans- 
ferred into the PAH in order to remove a carbon atom and repre- 
sents the analog of the displacement energy in a solid. For an 
extended discussion about the determination of the threshold en- 
ergy Tq we refer the reader to MJTb (Sect. 2.2. 1). We recall here 
that the value of To is, unfortunately, not well established, be- 
cause there are no experimental determinations on PAHs and the 
theoretical evaluation is uncertain. We decided to explore possi- 
ble values: 4.5 and 7.5 eV, close to the energy of the single and 
double C-bond respectively and 15 eV, compatible with the ex- 
pected threshold for a single walled nanotube. We adopt 7.5 eV 
as our reference value, consistent with all the experimental data. 

To calculate as a function of the kinetic energy of the inci- 
dent electron, it is important to remember that we are consider- 
ing relativistic particles, thus the appropriate expression for the 
kinetic energy is the following: 



1 



- 1 



(19) 



From Eq.[T9]we then derive fi 



nif. c- 



(20) 



The displacement cross section cr calculated from Eqs. [TSland 
I20I is shown in Fig.|4]as a function of the electron kinetic energy, 
for three different values of Tq. Above threshold the cross section 
increases with electron energy and decreases again at higher en- 
ergies because of relativistic effects, reaching the constant value 
given by the following asymptotic expression 



cr • 



%ZlE\nal 1 
M2c2 7^ 



(21) 



As expected cr decreases for increasing values of the thresh- 
old energy. Around the peak the change of To from 4.5 to 15 eV 
introduces a variation in the cross section of a factor of about 9, 
which reduces to 3.4 - the ratio 15/4.5, cf. Eq.|2T|- for electron 
energies above ~2 MeV. 

4. The CR spectrum 

The stopping power of ions with energy above ~1 MeV/nucleon 
decreases with increasing energy (cf. Sect. 2.1). This implies that 
the CRs responsible for the major energy transfer to a PAH, 
and subsequent damage, are the lower-energy ones (below 1 




Kinetic energy (IVleV) 

Fig. 4. High energy electron cross section for carbon atom re- 
moval, calculated for three values of the threshold energy To- 



GeV/nucleon). Unfortunately this part of the CR spectrum is 
not accessible fro m the Earth bec ause of the phenomenon called 
solar modulation (IShapirdll99lh . CRs entering the heliosphere 
see their intensity reduced by the effect of the solar wind, espe- 
cially at low energies and when the solar cycle is at its maxi- 
mum. The solar magnetic field is frozen within the plasma of the 
solar wind and drawn out with it into a spiral structure. CRs en- 
countering the solar wind are then convected outward. Moreover 
when such charged particles interact with the expanding mag- 
netic field, they are adiabatically decelerated. Hence, the CRs 
observed at a given energy were originally much more energetic. 

Because of the solar modulation, the interstellar CR spec- 
trum at low energies needs to be evaluated theoretically by solv- 
ing the transport equation for particles in the ISM, assuming 
an appropriate CR spectrum at the sources and taking into ac- 
count all possible mechanisms able to modify the intensity of 
the CRs during their propagation (energy losses, fragmentation 
etc. Shapiro 1991). The Pioneer and Voyager spacecraft have 
probed the heliosphere out to beyond 60 AU, greatly improving 
the understanding of the spectra of protons and heavier nuclei 
with energies above ~ 100 MeV and the effects of solar modu- 
lation, although limited information is available on CR nuclei 
below ~ 100 MeV d Webber. ,19981) . To describe the propaga- 
tion and escape of Galactic CRs (at energies below a few xlO'^ 
eV), a widely used approach is the leaky-box model, which as- 
sumes that the particles are confined to the Galaxy, with frequent 
visits to the disk boundaries where t hey have a smal l probabil- 
ity of leaking out (Ip & Axfordl9"85l : ISimpson & Garcia-Mun^ 



Il988l:llndrioloetal.ll2009h 

Following lWebber & YushakI (Il983h and lBringa et"an (l2007h 
the CR intensity I{E) as a function of the total energy per ion, E, 
is then given by 

/(£) = C£"-^^ / (e + E^f (cm^ s sr GeV)"' (22) 

The constant C can be determined by matching Eq. |22] with 
the high-energy CR spectrum measured on the Earth (see be- 
low). The scaling factor sets the level of low-energy CRs 
(IWebber & Yushaklll983h . which is ma inly determined by ion- 
ization loss and Coulomb collisions dip & AxfordI Il985h . At 
higher energies (above ~ 1 GeV/nucleon) diffusive losses domi- 
nate. 

The interstellar CR spectrum can be constrained by molec- 
ular observations. CR protons ionize atomic and molecular hy- 
drogen and this ionization drives interstellar chemistry through 
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Table 1. CR ion spectra parameters. 



Ion 






CM 


j(d) 

'o 




7 


H 


1.0 


0.12 


1.45 


11.5x10 


-9 


2.77 


He 


4.0 


0.48 


0.90 


7.19x10 


-9 


2.64 


CNO 


14. 


1.68 


0.36 


2.86x10 


-9 


2.67 


Fe-Co-Ni 


58. 


6.95 


0.24 


1.89x10 


-9 


2.60 



(a) 

(b) 

(0 
id) 



M\ in amu. 

£• (GeV) = £;(H)x Mi (ion). 
In units of (cm^ s sr GeV"''^)^' 
In units of (cm^ s sr GeV)"'. 



ion-molecule reactions. Analysis of molecular observations in 
diffuse clouds result in an average primary io nization rate of 
4 X 10"'^ s ' (H-nuclei) ' dlndriolo et al.ll2009l and references 
therein). The CR ionization rate follows from a convolution of 
the CR spectr um, Eq.[22] with the hydrogen ionization cross sec- 
tion dBringa et al..2007i) . The scahng factor E*q can be calculated 
then from 



^ = 5.85x 10""'(£(V0.1GeV) 



-2.56 



s ' (H nuclei) ' 



(23) 



Again, because of the steep decrease of the ionization cross sec- 
tion with energy, the CR ionization rate is most sensitive to the 
low energy CR flux. For hydrogen this results in iij ~ 0.12 
G eV. For heav i er par ticles we adopt the same scaling rule 
as iBringa etall (l2007h : El{\on) = £'*(H)x Mi GeV particle \ 
where M\ is the mass of the particle in amu. 

To calculate the constant C, we matched Eq. |22] with the 
high-energy spectrum detected on Earth. For the high energy 
data, which are not i nfluenced by the solar mod ulation, we adopt 
the expression from lWiebel-Sooth et al.l ( Il998l) 

/(£) = h [E (GeV)/(1000 GeV)]"^ (cm" s sr GeV)"' (24) 

where /q and y depend on the CR ion. 

In this study we consider the most abundant CR compo- 
nents: H, He, the group C, N, O and the group Fe, Co, Ni. 
The latter are in fact often detected as a group because of the 
experimental difficulty in distinguishing between particles with 
similar mass. The spectra were calculated using the method de- 
scribed above, with the h igh-energy parameters /q and y from 
IWiebel -Sooth et a The matching between the low and 

high energy regimes is at £ = 1 TeV A list of the parame- 
ters required for the calculation is reported in Table [T] and the 
resulting CR spectra are shown in Fig. |5] For the lowest en- 
ergy of the interstellar CR spectrum, we choose the value of 
5 MeV/nuc leon, consistent wit h the limit of ~ 1 MeV/nucleon 
adopted by lip & AxfordI (Il985h for their calculation of the CRs 
spectrum, and which corresponds to the lower limit of the energy 
range where ionization loss rapidly diminishes the propagation 
of CRs in the ISM. 

The same approach as that used for heavy particles (ions) 
can be applied to CR electrons. In this case solar modulation 
also alters the spectrum of the electrons entering the solar cav- 
ity. The interstellar spectrum at low energies has to be calculated 
solving the transport equation for electrons in the ISM, taking 
into account the energy loss processes relevant for electrons, i.e. 
bremsstrahlung, synchrotron and inverse Compton. The low en- 
ergy spectrum then has to be connected to the measured high- 
energy spectrum, which is not affected by the modulation. We 



adopt the expression from ICummingsl (1 19731). calculated in the 
framework of the leaky-box model (see also lip & Axfordll 19851 : 
Moskalenko & Strong 1998) 



I{E) = A [£(GeV) x 10^] (cm" s sr GeV)" 
where E is in GeV and 



(25) 



A = 0.0254 x 10\ r = 1.8 for 5 x 10"^ <E<2, 
A=5.19xl0^ 7 = 2.5 for2<£<10^ 

The calculated spectrum is shown in Fig. |5] For the lowest 
energy, we again adopt a value of 5 MeV, coherent with the range 
of influence of bremsstrahlung, the dominant energy loss mech- 
anism for electrons with energies less than a few hundreds of 
MeV. 

It is important to underline the fact that the method used 
to calculate the interstellar CR spectrum at low energies, i.e. 
solving the transport equation for particles adopting appropriate 
parameters in order to match the measurements, does not pro- 
vide an unique solution. It has been shown that different sets 
of parameters for the transport model result in a set of differ- 
ent spectral shapes all consistent with the CR intens ity mea- 
sured at higher energies (see e.g. iMewaldt et al.ll2004l) . Strong 
constraints on the low-energy CR spectru m have been provided 
by the interstellar i onization rate (e.g. Ilndriolo et al.l 120091 : 
iPadovani et alJl2009l and references therein), which all owed to 
exclu de some models, for instance the the W98 model ( Webbed 
119981) for protons and the COO model dMoskalenko et alJl2002l) 
for electrons. Despite such improvement, the ambiguity on the 
low-energy region of the CR spectrum has not been completely 
solved yet. This uncertainty does not affect the conclusions of 
our present work, but it is important to keep it in mind. 



5. Collision rate and C-atom ejection rate 

To calculate the collision rate between PAHs and CR ions and 
electrons, with energy fluxes described by Eq.l22landl25]respec- 
tively, we follow the procedure illustrated below. 

For ion collisions, where the energy is transferred to the 
whole PAH via electronic excitation, we first need to calculate 
the term 2(£'), which takes into account all possible ion trajec- 
tories through the PAH, with their corresponding transferred en- 
ergies and dissociation probabilities. Adopting the configuration 
shown in Fig.|2]we have 



Z(£) 



27r/ 



dQcr<,()5-) P{E,d) 



•nil 



o-M)P{E,&) sini9-di? 



(26) 



(27) 



i>=0 



with f2 = sin 1? di? d(p, with if running from to 2;7r and § from 
to n/2. The geometrical cross section seen by an incident particle 
with direction defined by the angle -& is given by 

o-g{§) = ttR^ cos§ + 2 Rd sm§ (28) 

which reduces to o-g - nR^ for )? = (face-on impact) and to 
(Tg — 2Rd for § — nl2 (edge-on impact). The term P{E, ff) rep- 
resents the total probability for dissociation upon ion collision, 
for a particle with energy E and incoming direction The ion 
collision rate is calculated by convolution of the term ^(ii) over 
the CR spectrum h{E) 



R\ 



,coll 
CR 



Fc/i(£)S(£)d£ 



(29) 
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Fig. 5. Interstellar CR spectrum of H, He, CNO, Fe-Co-Ni and 
electrons as a function of the particle energy. 



Because each (electronic) ion interaction leads to the removal of 
two carbon atoms from the PAH, to obtain the C-atom ejection 
rate the collision rate has to be multiplied by a factor of 2: 



CR 



Fc Ii(E) I.(E) dE 



(30) 



The factor, Fq, takes Coulombian effects into account (see be- 
low). For iimax we adopt a value of 10 GeV, corresponding to the 
highest energy for which experimental stopping determinations 
exist and thus for which Eq. |2]is valid. Moreover the CR inten- 
sity and stopping power decrease rapidly, so ions with energy 
above 10 GeV do not contribute significantly to the integral in 
Eq. [29] Concerning the lower integration limit, we perform the 
calculation for = 5 MeV/nucleon, the lower limit assumed 
for the CR spectra. 

For interactions with CR electrons, each binary collision re- 
sults in the ejection of one single C-atom. Thus the ejection rate 
coincides with the collision rate, and are both given by the fol- 
lowing relation 



e, CR 



R 



e,CR 



0.5A?c47r 



r 



Fch(E)cr(E)dE 



(31) 



where the ejection cross section per target atom cr has to be mul- 
tiplied by the number of C-atom in the molecule, A^c- The factor 
0.5 takes the angle averaged orientation into account (MJTa). As 
for the CR ions we assume as integration limits the values E„iin 
= 5 MeV and ^max 10 GeV. For the electrons, the upper limit 
is not constrained by the stopping theory (the expression for the 
ejection cross section holds for even higher energies) and is only 



related to the steepness of the spectrum which results in a negli- 
gible contribution from electrons with energy above 10 GeV (cf. 
Fig.EJ. 

The Coulombian correction factor Fc (MJTb) takes into ac- 
count the fact that both target and projectiles are charged, and 
that the collision cross section could be enhanced or diminished 
depending on the charge of the PAH. Because we are consider- 
ing high energy ions and electrons, which are unaffected by the 
Coulombian field, Fc is always unity. 

6. Results 

6.1. PAH lifetime 

Collisions with CR ions and electrons will cause a progressive 
decrease in the number of carbon atoms in a PAH molecule. For 
interactions with ions, after a time t this number is reduced to 



A^c(0 = A^c(O) - RicRt 

and the fraction of carbon atoms ejected from this PAH is 

^i,CR t 



A^c(O) 



(32) 



(33) 



where Ri cr is the C-atom ejection rate from Eq.[30] We assume 
that the PAH is destroyed after the ejection of 1/3 of the carbon 
atoms initially present in the molecule. This occurs after a time 
To which we adopt as the PAH lifetime against CR ions bom- 
bardment, and is given by 



To = 



Nc 
3Ri, 



CR 



(34) 



For electron collisions, the number of carbon atoms in the PAH 
molecule after a time t is 

A?c(0 = A^c(0)exp[-r/T'] (35) 

and the fraction of carbon atoms ejected from this PAH is 

FL(f) = (l-exp[-r/r']) (36) 

with the time constant t' = A^c /Re, cr, where 7?e, cr is the ejection 
rate from Eq.|3T] The ejection of 1/3 of carbon atoms originally 
present in the PAH molecule, after which the PAH is considered 
destroyed, takes the time to given by 



(37) 



As for ion collisions, we adopt to as the PAH survival time 
against CR electrons. The ejection rate Re,cR scales linearly with 
A^c, hence the corresponding lifetime to is independent of PAH 
size. 

For collisions with ions, we calculated the lifetime of PAHs 
of four sizes Nc = 50, 100, 200 and 1000 (radius apah = 6.36, 9, 
12.7, 28.5 A, «niax = 10, 20, 40 and 200 respectively) bombarded 
by CR H, He, CNO and Fe-Co-Ni ions, assuming three different 
values for the electronic binding energy, £0 = 3.65, 4.58 and 
5.6 eV (cf. Sect. 4.1 in MJTb). In principle, not all the energy 
transferred to the PAH will be internally converted into vibra- 
tional modes, with consequent relaxation through dissociation 
(or IR emission). Other processes can occur, for instance the 
production of Auger electrons which will carry away from the 
molecule a part of the transferred energy. We don't know ex- 
actly how quantify this energy partitioning and so we introduce 




Fig. 6. PAH survival time against collisions with CR ions, as a function of the fraction / of the transferred energy that goes into 
vibrational excitation. The lifetime has been calculated for PAHs of different sizes (Nc), assuming three values for the binding 
energy of the ejected fragment, Eq. The line style indicates the type of ion: dotted for H, dashed for He, solid for CNO and dot-dot- 
dashed for Fe-Co-Ni. The line color indicates the value of the fragment binding energy; gray for £0 = 3.65 eV, black for Eq - 4.58 
eV and red for £0 = 5.6 eV. 



a factor /, which represents the fraction of the transferred energy 
T that goes into vibrational excitation. 

The survival time as a function of the factor / is shown in 
Fig. |6] The PAH lifetime becomes shorter as more energy goes 
into vibrational excitation (increasing /) and for lower values of 
£0, because this implies a higher temperature, T-^y, and a lower 
energy to eject a fragment, resulting in a larger dissociation prob- 
abiUty (Eq.inii. 

The dissociation probability is more sensitive to £0 when 
the energy available for dissociation is lower, i.e. light projec- 
tile and small /, and when the same amount of energy has to 
be spread over more bonds (increasing size). This explains why 
the separation between the time constant curves corresponding 
to the different values of £0 decreases with increasing available 
energy and mass of the projectile (from H to Fe-Co-Ni), while 
it gets bigger for larger PAHs. Big PAHs are more resistent to 
CR bombardment because for any given transferred energy their 
dissociation probability is lower. 



The lifetime against Fe-Co-Ni bombardment is essentially 
constant - except for very low available energy and very large 
molecules (A^c = 1000). Its large value, tq ~ few 10^ yr, results 
from the fact that, despite the huge amount of energy transferred 
into the molecule (cf. Fig. [T]i, the Fe-Co-Ni abundance in CRs 
is small. This leads to a low collision rate and long lifetime. 
For hydrogen, the high abundance is not enough to compensate 
for the small stopping power, which rapidly decreases above 1 
MeV/nucleon. As a result the lifetimes are long. Since the en- 
ergy transferred to the molecule is small, the collision rate and 
the survival times are sensitive to the adopted parameters: the 
fraction of the energy transferred available for dissociation /, 
the fragment binding energy £0 and the PAH size. This variabil- 
ity is shown in Fig. |6] The situation for helium and CNO lies 
somewhere between the cases for H and Fe-Co-Ni. 

The CR electron time constant is independent of the PAH 
size, and has been calculated for three values of the threshold 
energy for carbon atom ejection. To = 4.5, 7.5 and 15 eV. The 
results are shown in Fig. |7] and summarized in Table |2] The 
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Table 2. Time constant tq for carbon atom ejection following 
PAH collisions with CR electrons. 





Ta (eV) 






4.5 


7.5 


15. 


1 £™„(MeV) 


To(yr) 1.2x10'^^ 


2.0x10'^ 


4.2x10''' 


1 5 



PAH lifetime is longer for increasing values of Tq because of 
the diminution in the cross section cr. Given that cr is almost 
constant for energies above ~ 2 MeV, for a fixed To this part of 
the spectrum does not contribute to the variation. Because of the 
small ejection cross section and the steep CR electron spectrum, 
all the calculated time constants are long (tq >10''' yr). 

6.2. Discussion of the uncertainties 

The first source of uncertainty affecting the study of PAH inter- 
action with high energy particles is the accuracy of the stopping 
theory used to calculate the energy loss of ions in matter. This is 
a difficult issue because the current theory, based on the Bethe- 
Bloch equation and described in Sect. 2.1, is in fact a combina- 
tion of different theoretical approaches with corrections coming 
from fits to the experimental data. A better way to pose the prob- 
lem is by instead asking: "How accurately can stopping pow- 
ers be calculated? " The comparison between theory and experi- 
ments, together with the evaluation of possible variation sources 
such as structural variations in the targets, provides an accuracy 
ranging from 5 to 10 % (Ziegler 1999). 

In our calculation we do not use the complete Bethe-Bloch 
equation (implemented in the SRIM code with all corrective 
terms) but the analytical approximation given by Eq. [10] The 
discrepancy between the stopping power curves from these two 
formulations is very small, and introduces an uncertainty of at 
most 10 % in our calculation. An additional source of uncer- 
tainty comes from the calculation of the energy transferred to 
the PAH using the approximation in Eq. [TT] which assumes a 
constant stopping power along the distance travelled through the 
molecule. In this case the difference with the exact calculation is 
also limited to less than 10 %. 

Concerning the collisions with high energy ions, the domi- 
nant source of uncertainty in the determination of the PAH life- 
time is the fragment binding energy £0 (cf. MJTb, Sect. 4.1), 
which is strongly modulated by other parameters: type of ion, 
PAH size and fraction of energy available for dissociation, as 
clearly shown in Fig.|6] To give an example, when Eq goes from 
3.65 to 5.6 eV, the time constant for a 50 carbon atoms PAH col- 
liding with helium varies by a factor of ~ 2x10' for / = 0.25 and 
by a factor of 8 for / = 1 . For CNO impacting on a 200 C-atom 
molecule the change is a factor of 7 and 3 for / - 0.25 and / 
= 1 respectively. These numbers give an idea of the huge and 
complex variability induced by the uncertainty in the parameter 
Eo. 

In the treatment of PAH collisions with high energy elec- 
trons, two sources of uncertainty have to be considered: the an- 
alytical approximation to the numerical solution for the ejec- 
tion cross section (Eq. [TSl l and the choice of the threshold en- 
ergy T(). While the discrepancy between the analytical expres- 
sion and the numerical solution for cr is only a few percent 
(McKinlev & Feshbach 1948i), the variation of Tq from 4.5 to 
15 eV implies a variation in the cross section of about a factor of 
9 around the peak (electron energy ~ 0. 1 MeV), and a factor of 3 
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Fig. 7. Time constant tq for carbon atoms removal due to col- 
lisions with CR electrons, as a function of the threshold energy 
7^0. 



above 2 MeV, where the cross section becomes almost constant. 
This results in a change of the electron time constant of a factor 
of ~ 3.4, consistent with the cross section variation. To summa- 
rize, in our study we explored a wide range in values for To and 
Eq, and the resulting variation is in some sense the maximum 
possible. 

A final remark concerns the CR spectrum at low energies. As 
previously mentioned, the major contribution to PAH destruction 
comes from low-energy CRs, both because the energy transfer is 
more efficient for low energy projectiles (ions) and the CR in- 
tensity decreases with increasing energy (for both ions and elec- 
trons). Unfortunately, because of the solar modulation this part 
of the spectrum is not accessible to measurements from Earth 
and needs to be evaluated theoretically. The result of this evalua- 
tion is not unique, and various spectral shapes of the interstellar 
CR spectrum have been shown to be consistent with the mea- 
sured CR intensity (cf. Sect. 4). The spectra we adopt, which 
of course match the measured high energy spectra, have been 
calculated using a rigorous physics and taking into account all 
possible constraints from available measurements, including the 
stringent requirements from the CR ionization rate. Nevertheless 
the lack of data in such a crucial CR energy region still repre- 
sents a source of uncertainty that is important to keep in mind, 
although it does not have a significant impact at the level of ac- 
curacy for our study. For an overview of the current state of the 
art of the uncertaintie s in the CR fluxes we refer the reader to 
iPadovani et al.l (l2009l) and references therein. 

7. Discussion 

7. 1 . The lifetime of PAHs against CR processing 

In Sect. 6.1, we estimated the PAH lifetime against CR bom- 
bardment in the neighborhood of the Sun, i.e. considering an in- 
terstellar spectrum based on galactic CR measurements near the 
Earth but corrected for the influence of the Heliosphere (cf. Sect. 
4). Fig.[8]shows the PAH lifetime against CR bombardment (ions 
+ electrons) as a function of the size of the molecule, Nc, for / 
- 0.5 and 1, assuming our reference values of 7.5 and 4.6 eV for 
the parameters To and £0 respectively. These lifetimes should 
be compared to the lifetime against shock destruction of PAHs 
of approximately 150 Myr (MJTa). Depending on the amount 
of transferred energy available for dissociation (50 or 100 %), 
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Fig. 8. PAH survival time against CR bombardment (ions + elec- 
trons) as a function of the molecule size (A^c)- The total lifetime 
has been calculated for / = 0.5 and 1, where / is the fraction 
of the transferred energy available for dissociation, and adopting 
our reference values for the threshold energy for carbon atom 
ejection, To = 7.5 eV and for the fragment binding energy, £0 = 
4.6 eV. We remind the reader of the variation in the calculated 
survival time against CRs, due to the uncertainty on the param- 
eters £0 and To (cf. Fig.|6]and Fig.|7ll. The PAH lifetime against 
shock destruction in the ISM is shown for comparison. 



CRs are calculated to be the dominant destruction agent in the 
warm ISM for PAHs with less then -150 and -280 C-atoms 
respectively. For larger PAHs, shocks take over. In this assess- 
ment, it should be kept in mind that shock processing occurs 
predominantly in warm intercloud medium of the ISM: shocks 
faster than 100 km s ' are very rare in the cloud phase since 
supernova remnants predominantly expand in the warm or hot 
intercloud medium and the shock speed in clouds is then down 
by the square roo t of the density ra tio between these different 
phases in the ISM (I Jones et al ][T99l . Shock processing of PAHs 
in diffuse clouds will then predominantly occur when cloud ma- 
terial exchanges with the warm intercloud medium (where shock 
timescales are short). The timescale for this exchange is included 
in the lifetime estimate for shock processing. CRs on the other 
hand penetrate all phases of the ISM, except perhaps the dens- 
est molecular cloud cores, and can process all material that they 
interact with. 

Finally, we note that the PAH survival times against CR elec- 
trons (to >10'^ yr) are longer than the Hubble time, implying 
that CR electron collisions are not important for the processing 
of interstellar PAHs. We point out the fact that this study does 
not include the secondary electrons with an energy of about 30 
eV produced by primary CR electrons by the initial ionization of 
H or H2 in diffuse clouds. Such electrons are energetic enough 
to induce some PAH fragmentation but, although initially abun- 
dant, they decay easily. 



7.2. PAHs in galactic Inaios 

ISO and Spitzer observations revealed that PAHs are abundant 
in the halos of normal spiral galaxie s (I Irwin & Maddenll2006t 
llrwin et al 1 120071: IWhalev et al.ll2009l) . Here, we will discuss the 
potential role of CRs in the destruction of these PAHs. 

After having estimated the PAH lifetime against CR bom- 
bardment in the vicinity of the Sun, the question that needs to be 
addressed is whether the CR spectrum changes across the galaxy 



and with the galactic latitude, and how this may affect the PAH 
lifetime. 



7.2.1. The CR spectrum in galactic lialos 

A powerful tool to probe the presence of CRs in regions distant 
from the Solar System has been provided by gamma and radio 
astronomy. CR protons interact hadronically with the nuclei of 
the interstellar medium (mainly pro tons and a particles) produc - 
ing both charged and neutral pions (iFermill 1 950t ISteckedl 1 973h . 
Charged pions decay into secondary leptons, while neutral pions 
decay into two y photons: 



p+N 

7T~ 



2y 



(38) 
(39) 
(40) 



Electrons loose energy via synchrotron and inverse Compton, 
producing detectable radio emission, while y-rays are directly 
measurable by gamma-ray telescopes. The quantity of cosmic- 
ray protons can then be calculated from their decay products, 
moreover, because photons are not affected by magnetic fields, 
they can pro vide a direct indication about the location of their 
sources (e.g. iRamana Murthv & Wolfendaldl 1 9931: [Hunter et all 
[T997I) . 

Different studies (see IShibata et al.l l2007l and references 
therein) have shown that the slope of the energy spectrum is al- 
most independent of the observational site in the Galaxy, but this 
is not true for the amplitude of the spectrum. In the galactic plane 
the C R intensity shows a longitudinal gradient (IShibata et al.l 
l2007h . with an exponential decrease from the Galactic center 
toward the periphery with a scale length of 20 k pc. Assuming 
the galactic parameters from IShibata et al.1 (l2007h we obtain an 
enhancement at the Galactic center of about a factor 1.5 with 
respect to the Solar System. 

Radio and gamma observations, t ogether with the study 
of unstable CR nuclei such as '"Be (iHavakawa et al.l 119581: 
ISimpson & Garcia-Munozlll988l) have demonstrated that CRs 
are not confined to the thin disk w hich contains their sources, 
as ass umed by the leaky-box model (ISimpson & Garcia-Munozl 
Il988l) . but are able to travel long distances outside the 
galactic plane . Various models for CR propagation (e.g . 



Ginzburg et all Il980t iGinzburd 119881: feerezinskii et"an 1 19901: 
Strong & Moskalenkol 19981) assume a thin disk (half-thickness 



hi ~ 100 - 200 pc) located on the galactic plane, where the 
sources reside, surrounded by a large diffusive halo with half- 
thickness /zh ~ 1 - 15 kpc , where CRs spend part of their lif e. 

Using the formalism developed by IShibata & ltd (|2007|) the 
latitudinal gradient of the CR intensity for a given ion / can be 
written as 



Ii(r,z,E) 



^ exp(-|z|/zD) 



(41) 



where r is the radial distance from the center of the galaxy pro- 
jected on the galactic plane and z is the latitudinal distance from 
the Galactic plane. Eq.|4T]shows that the CR intensity attenuates 
with the scale height zd of the diffusion coefficient. The lati- 
tudinal gradient is calculated assumi ng the value zd = 2.5 kpc 
proposed by e.g. IShibata et al.l ( l2007h and is shown in Fig.|9] 

Ultimately, CRs are thought to be accelerated by supernova 
remnants, tapping some 10% of the SN energy. For other galax- 
ies, we will scale the overall CR density with the star forma- 
tion rate of the galaxy and we will adapt scale lengths and scale 
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Fig. 9. The latitudinal CR gradient in both Milky Way (MW) and 
NGC 891, as a function of the vertical distance z from the galac- 
tic plane together with the vertical variation of the interstellar 
radiation field Go in NGC 891 [courtesy F. Galliano). The CR 
intensity / decreases with the scale height zd of the diffusion 
coefficient. 



heights appropriate for the Milky Way. Since the CR intensity 
decreases exponentially outside the galactic plane, we expect the 
R\H lifetime at higher latitudes to be enhanced by the same fac- 
tor. 



7.2.2. The origin and evolution of PAHs in galactic halos 

PAH emission has been observed at high la titudes in several 
normal galaxies with scale heights of 2-3 kpc (I Irwin & Madded 
[2006; Irwin et al. 2007; Whaley et al. 2009J. The study of NGC 
891, the edge-on twin of the Milky Way, is particularly instruc- 
tive. While the vertical distribution of the 8 fim PAHs emission 
is noticeably narrower than that of the cold dust radiating at 450 
jum, modeling reveals that this likely reflects the variation of 
the stellar population from the disk to the halo; e.g., the dom- 
inant (FUV) heating sources of the PAHs are located in the disk 
while the old stellar population, which contributes substantially 
to the heating of th e dust, extends some 3 kpc above the disk 
(IWhalev et al.ll2009h . 

Some of the PAH molecules, as well as the dust, present at 
high galactic latitudes may originate from mass losing AGB stars 
in the halo. For the Milky Way, recent optical and infrared sur- 
veys (2MASS and DENIS) have revealed the presence of C-rich 
AGB stars ensh rouded in dusty ejecta at distances of 2-6 kpc 
above the plane dGroenewegen et al.ll997l:[Mauronk2008.) . Some 
100 such C-stars are known in the halo up to distances of 150 
kpc and man y of these belong to the tidal stream of the Sgr 
dwarf galaxy (iMauron & Kendalill2005l) . There is no complete 
census of such objects and a mass balance of PAHs and dust in 
the halo cannot yet be assessed. 

More likely, though, most of the PAHs and dust represent 
signpost of the large scale circulation of matter between the disk 
and halo. Indeed, the vertical distribution of PAHs is similar to 
that of the dust - once the differences in heating are accounted 
for - and to other tracers of this circulation pattern such as the 
diffuse ionized gas. If PAHs are indeed transported to high lati- 
tudes through the action of the galactic fountain, then the accel- 
eration associated with this process must have been very gentle 
since PAHs are quickly destroyed in shocks faster than 100 km 
s (MJTa). Likely, PAHs are entrained as cloud(let)s sheared 



off the chimney walls in the hot gas of the venting supemovae 
and lifted high above the plane. As shown in MJTb, PAHs are 
rapidly destroyed in hot gas. Hence, after evaporation of these 
cloudlets driven by thermal conduction from the hot gas, PAHs 
will be gone and clouds subsequently condensing through ther- 
mal instabilities in the halo will not show up in the PAH emission 
bands. 

Hence, PAHs in the halo may represent a 'molecular dye' 
with which this entrainment aspect of the disk-halo circulation 
flow can be separated from the evaporation/condensation pattern 
and followed in detail. The effectiveness of this PAH-dye is mod- 
ified on the one hand by the diffusion of the FUV photons needed 
for their excitation from the disk and on the other hand by de- 
struction of PAHs in the harsh environment of the halo. Here, we 
are concerned with the latter aspect, the lifetime of PAHs at high 
latitudes above the plane. First, the residence time is of inter- 
est. With a total mass of Warm Intercloud material of ~10'' Mq 
in the lower (~ 0.5 kpc) halo and a circulation rate of 5 M0/yr 
between the plane and the halo for the Milky Way, the residence 
time of PAHs at these latitudes is some 2x10*^ yrs. The residence 
time at the higher latitudes (~ 2 kpc) we are interested in here 
will be commensurately larger. At the latitudes where PAHs are 
observed (~2 kpc), we expect that supernova shock waves, an 
important destruction agent in the disk of the galaxy, are of lit- 
tle concern. UV photolysis is generally considered to be a main 
agent for the destruction of small (<50 C-atoms) PAHs, weeding 
out the less stabl e (e.g., smalle st and/or non-compact) PAHs on a 
rapid timescale (lTielensll2008l) . Indeed, the minimum size in the 
PAH-size-distribution is thought to reflect this process. However, 
because this process is very sensitive to size, compact PAHs only 
slightly larger than this minimum size are essential 'indestruc- 
tible'. Thus, we expect that PAHs transported upwards from the 
plane by the galactic fountain have already been weeded down to 
the most stable forms and because of the decreased UV flux with 
latitude (Fig.|9]l, further UV photolysis will be of little concern. 
Using the results of this paper, we can estimate the destruction 
of PAHs by CRs. 

Taking the CR distribution in the halo from Sect. 7.2.1 we 
calculate the CR destruction timescale. Because of its similarity 
with the Milky Way, we adopt for NGC 891 the same diffusion 
scale height zd - 2.5 kpc. This implies that the latitudinal gradi- 
ent from Eq.|4T]is the same for both galaxies. The CR intensity 
variation in NGC 891 with respect to the solar neighborhood is 
given by the following expression 



W89l(0,k|,£) ^ 3.8^ , , ^ 

— — — = (1.5 X — ) exp(-|z|/zD) 

4mw('"o, 0, E) 3 



(42) 



i 



where 1.5 is the CR enhancement factor in the galactic center 
with respect to the solar neighborhood and the factor 3.8/3 takes 
into account the increased cosmic-ray intensity in NGC 891 due 
to its higher star formation rate (3.8 Mq yr jPopescu et al 
2004 with respect to the Milky Way (~3 M© yr^MScalo 198 _ 
Prantzos & Silidl998h . At the typical PAH detection scale height 
|z| - 3 kpc we obtain a decrease of the CR intensity of a factor of 
0.6. The lifetime against destruction will be then enhanced by a 
factor of (1/0.6). 

The total PAH lifetime in the halo of NGC 891 compared 
with the circulation timescale is shown in Fig. \W\ Depending 
on the amount of energy available for dissociation, the survival 
time against CRs for PAHs with more than ~100 C - atoms is 
comparable or longer than the circulation lifetime. PAHs could 
be good tracers of material entrained in the plane and carried 
to high latitudes by the galactic fountain as long as this mate- 
rial never loses its 'identity' through evaporation into the hot gas 
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Fig. 10. Total PAH lifetime against CR bombardment in the halo 
of NGC 891, evaluated at an altitude of 3 kpc from the galactic 
plane, compared with the circulation timescale between the disk 
and the halo. The value of the circulation lifetime represents a 
lower limit, as indicated by the arrow. Otherwise as Fig.|8] 



followe d by recondensation i nto cloudlets. Of course, as empha- 
sized bv lWhalev et"an (l2009l) . such a study has to properly eval- 
uate the flux of the pumping of UV photons in order to be of 
quantitative value. In that case, mid-lR studies can potentially 
trace this entrainment process and the exchange between cool 
and hot phases on arcsecond size scales. Nevertheless other fac- 
tors have to be taken into account, such as the possibility of local 
PAH production in the halo. 

7.3. PAHs in galactic winds 

Starburst galaxies are characterized by very intense star forma- 
tion, usually concentrated in the nuclear region of the galaxy, 
which can occur as a single instantaneous burst or as a high star 
formation activity sustained over a period of time. Since galactic 
supernovae are considered the most probable acceleration source 
for CRs wit h energy below the "knee" (few 10'^ eV - see e.g. 
iDrurvl [19941) . we expect that an enhancement of the SFR will 
result in a corresponding increase of the CR intensity and dust 
destruction by shocks. 

In our previous work (MJTb) we discussed the observation 
of PA H molecules in the bipolar outflow of th e starburst galaxy 
M82 (lArmus et al ] l2007HGalliano et al.ll2008h . PAHs have been 
detected far outside the galactic p lane at least up to 6 kpc from 
the disk ( lEngelbracht et al.||2006b . in a region which is spatially 
correlated with an extended X-ray emission originating from 
tenuous hot gas (electron density «e = 0.013 cm"^, temperature 
T = 5.8x10^ K). Our study has demonstrated that, under these 
conditions, PAHs are destroyed by collisions with thermal elec- 
trons on a timescale of few thousand years. Our conclusion was 
that the PAH survival is possible only if the molecules are iso- 
lated from this hot gas, probably in cooler cloudlets of PDR-type 
gas entrained in the galactic wind. 

In fact supernova-driven galactic-scale winds, often called 
"superwinds", are comm only observed in starburst ga laxies, 
including M82 (see e.g. [Strickland & HeckmanI |2009|. and 
references therein). In addition, theoretical studies jlpavichl 



1975 ^. 'Breitschwerdt et al.l ll99lL 11993 1: 'Zirakashvili et al 1996^ 
Ptuskin et al. 1997), together with radio (Heesen et al. 2009) and 



port of CRs outwards into the halo, in addition to the diffusive 
transport predicted by CR propagation models (cf. Sect. 7. 1). 

Assuming also that the M82 superwind could be partially 
CR-driven, we can try to evaluate the CR intensity at the high 
latitudes where PAHs are detected, to assess PAH lifetimes 
against CR bombardment and compare them with the survival 
time in hot gas of the outflow. Unfortunately the properties of 
superwinds at high galactic latitudes are not well constrained. 
Theoretical studies mainly concern small regions around the nu- 
cleus of the galaxy (radii of few hundred pc from the ce nter), 
where the transport of matter is dom inated by advection jVoUd 
20031: IStrickland & Heckmanll2009h . while in the outer regions 
diffus ion is expected to be dominant (e.g. IZirakashvili & Volkl 
l2006l for the starburst galaxy NGC 253). This is not the case in 
the Milky Way, where the CR transport is diffusive on a kpc scale 
and b ecomes convective only at larger distances dPtuskin et al.l 

On the other hand, the CR propagation in M82 is not yet 
totally understood. Star formation, and thus SNRs which re- 
presents the main CR sources, are concentrated in the small nu- 
clear region from which the galactic superwind originates, we 
could then expect CRs being advectically transported by the 
wind out of the galactic plane up to a certain height (~ few hun- 
dred pc), and diffusing after in the volume occupied by the su- 
perwind, which is roughly a cylinder flaring at high z with height 
of ~ 12 kpc on the North side and 7.5 kpc on the South side 
( lLehnertetal.il 19991: IStevens et aT]|2003l) . It would also be pos- 
sible that CRs occupy a spheroidal halo as in the Milky Way, 
maybe diffusing back to the disk after leaving the wind region. 

Because of this lack of information, especially on the con- 
finement volume for the CRs, what we present here is a toy- 
model based on some simplifying assumptions, but which could 
be useful to interpret PAH observations in the outflow of M82. 
To calculate the CR spectrum in M82 we adopt the same spec- 
tral index as in the Milky Way (y = 2.7) and we assume that in 
the starburst nucleus the CR intensity scales with the supernova 
rate vsn (SN yr"') which is directly related to the star formation 
rate in the galaxy. Taking the values vsn(M82) - 0. 11 SN yr~' 
(iHuang et alJfT99l and vsn(MW) = 0.02 SN yr ' (iDiehl et al.1 
|2006|) we obtain / Imw(E) = vsN(M82)/ysN(MW) = 

0.11/0.02 = 5.5. 

Since PAHs are detected at vertical distances which are large 
compared with the typical height of the convective zone, we ne- 
glect this latter and assume pure diffusive CR transport above the 
galactic plane. We adopt the exponentially-decr easing CR latitu- 
dinal gradient from Eg. |4T] (tShibata et alJl2007l) with zd = /Jh/3, 
where /jj, is the maximum distance that CRs can reach above the 
plane through diffusion. Since the exact confinement region in 
M82 is not known, we just take equal to the vertical height 
of the galactic superwind (12 kpc in the North side, 7.5 in the 
South side), i.e. we assume that CRs can occupy at least as the 
same volume than the gas from the outflow. 

The PAH lifetime against CR bombardment at the vertical 
distance |z| above the starburst nucleus of M82 is then given by 
the following expression 



T0,M82(0, \Z\) = 



1 



To.MwjrQ, 0) 

5.5 X 1.5 exp(-|z|/zD) 



(43) 



X-ray observations (lEverett et alj|2008l) suggest the possibility 
of a CR-driven galactic wind, i.e. the possibility of "bulk" trans- 



where 5.5 and 1.5 are the two enhancement factor of the CR 
intensity in M82 with respect to the Milky Way (5.5) and in 
the galactic center with respect to the solar neighborhood (1.5), 
which will imply a decrease of the same factor for the PAH life- 
time with respect to the survival time ro,Mw('"o, 0) in the vicinity 
of the Solar System (Fig. [8]l. On the other hand the latitudinal 
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Fig. 11. Total PAH lifetime against CR bombardment in the out- 
flow of the starburst galaxy M82, evaluated at an altitude of 6 
kpc from the galactic plane. The PAH survival time is compared 
with the circulation, starburst and X-rays timescales. The values 
of the circulation and X-rays timescales represent lower limits, 
as indicated by the arrows. Otherwise as Fig. [8] 



gradient will result in a corresponding increase of to,m82 for in- 
creasing |z|, with scale height ZD(North) = 4 kpc and ZD(South) 
- 2.5 kpc. 

Let us calculate the PAH lifetime at the maximum distance 
|z| - 6 kpc where the molecules have been detected. This is done 
by scaling the time constants ro_Mw('"o, 0) in Fig. [8]by the fac- 
tor A ~1.3 (averaged over the two hemispheres) and the result 
is shown in Fig. [TT] The PAH lifetime against CRs is longer 
than the starburst timescale (tens of Myr), except for the smallest 
molecules when f - I, and for large PAHs (A^c > 150-200) it is 
comparable or longer than the circulation lifetime (~ 200 Myr). 
In MJTb we showed that in the outflow of M82 PAH destruc- 
tion due to X-ray absorption would take place on a timescale 
of at least ~ 20 Myr (assuming that each X-ray photon absorp- 
tion will lead to fragmentation). This is comparable to or larger 
than the starburst lifetime. PAHs dispersed in the hot gas will be 
rapidly destroyed by collisions with electrons. To survive they 
need to be isolated from this hot gas, presumably in a cooler and 
denser gas entrained in the superwind. On the other hand CRs 
and X-rays can penetrate in those clumps, thus the PAH lifetime 
there will be determined by the cosmic radiation and/or X-rays, 
depending on the PAH size and on the effective value of the sur- 
vival time against X-rays (cf. § 6.1 in MJTb). 

7.4. PAHs in galaxy clusters 

We consider now the PAH processing by CRs in an extragalac- 
tic environment, focusing on the specific case of cooling flow 
galaxy clusters. The measurements of high-energy gamma and 
neutrino radiation produced by the interaction of CRs (mainly 
protons) with the intra-cluster medium (ICM) have demonstrated 
that galaxy cluster are able to confine CRs for cosmological 
times (|Berezinskv et al. 1997). Our choice is thus motivated by 
the reasonable expectation to find a non-negligible population 
of CRs in such environments, although the spectral index and 
amplitude are not well constrained by observations. 

Many clusters of galaxies with cooling flows host a central 
galaxy which often shows nuclear activity. Many papers support 
the presence of CR protons in the relativistic plasma bubbles pro- 
duced by Active Galactic Nuclei (AGN), where the energetics 



are s ufficient to inject a significa nt amount of CRs into the ICM 
(see iPfrommer & EnBlinI |2004 and references therein). The 
slope of the injected spectrum depends on the fraction of par- 
ticles released from the plasma. If all the CR content is ejected 
by the radio plasma, we expect to find a moderately flat spectrum 
(ainj ~ 2.5) consistent with the flat spectra of CR electrons in- 
dicated by the radio emission from radio galaxies. On the other 
hand, if only a small fraction of the CR population is able to 
leave the plasma via diffusion, we expect to find an even flatter 
spectrum (ainj ~ 2.2) because the escape probability increases 
with momentum ( lEnBhnll2003h . 

Another possible source of CRs in galaxy clusters is repre- 
sented by the galactic winds from the central galaxy transport- 
ing CRs accelerated within the galaxy itself. Because Supernova 
Remnants (SNRs) are believed to be the ma in sources of 
Galactic CRs (e.g. lGinzbu"r3ll993l: lDrurvlll994l) . the expected 
injection spectral index is ainj ~ 2.4, if no further re-acceleration 
occurs. Apart from protons, the principal constituents of Galactic 
CRs are helium, carbon and iron ions. 

The CR population in galaxy clusters will then result from 
the combination of the above-mentioned sources, with a spec- 
tral index which reflects its composite origin. After injec- 
tion from the source, CRs start diffusing away from the cen- 
tral AGN. For the treatment of the transport of CR protons 
through the thermal ICM, we refer to the model elaborated by 
Pfrommer & EnBlin (2003)- The CR proton distribution function 
is described by a power-law in energy (or momentum) where the 
injection spectral index is modified by the diff'usion through the 
intra-cluster gas 

(r E) = f ^ V' (^y"" (44) 

GeV (/i-o'kpcj IceV/ 

with ffp = ainj + ffdiif ■ The coefficient aditf results from the fact 
that diffusion is momentum-dependent, and equals a; 1/3 for ac- 
tive diffusive transport of CR protons i n a Kolmogorov-like spec- 
trum of hydro-magnetic turbulence ( lKolmogorovlll94l1) . The 
model assumes the standard ACDM cosmology with Ho - 70 
li-io km s"' Mpc"', where /170 indicates the scaling with Hq. The 
normalization factor hcr.o is the crucial parameter required to 
det ermine the CR spectrum. 

IPfrommer & EnBlinI (|2004 constrained «cr,o for a sample 
of nearby cooling flow galaxy clusters with the aid of y-ray ob- 
servations, using the fact that CR protons interact hadronically 
with the thermal intra-cluster gas producing both charged and 
neutral pions (see Sect. 7.2.1). Charged pions decay into sec- 
ondary electrons/positrons + neutrinos/antineutrinos, while neu- 
tral pions decay into two y's [jfi 2y). The amount of cosmic 
protons can then be calculated from their decay products, which 
are detectable in the y-ray and radio bands. 

It is important to remember here that the threshold energy 
for pion production is Zsth = 0.78 GeV, implying that only 
the proton population with energy exceeding the threshold is 
able to produce pions hadronically, and can then be constrained 
by y-ray and radio observations. At lower energies the domi- 
nant energy loss mechanism for CR protons is Coulombian dif- 
fusion on the electrons in the plasma (electronic excitation), 
which results in the depopulation of the CR energy distribu- 
tion with consequent modificatio n of the power-law spectrum 
(iMannheim & Schlickeiseijll994 . 

Unfortunately, it is not possible to observationally constrain 
the low-energy CR pop ulation in galaxy c l usters, thus follow- 
ing the argument from iNath & BiermannI (Il994h we decided 
to extrapolate Eq. |44] down to a minimum energy E^i^ rep- 
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resenting the characteristic energy scale at which the spec- 
tmm deviates from a power-low. A similar approach is used by 
[Pfrommer& Enfilin (2004) as well. For the cut-off energy we 
adopt the value E ^m - 50 MeV, con sistent with the range ~ 50 - 
75 MeV found bv lNath et all (|2006|) in their study of the produc- 
tion of ^Li by CR protons in galaxy clusters. iSmin represents the 
lower energy that an energetic proton has to reach before suffer- 
ing large Coulombian losses, and can be equivalently expressed 
in terms of the quantity of matter to be traversed in order to ex- 
perience the same Coulombian losses. 

As previously mentioned, we assume the presence in the 
cluster CR population of the main components of the Galactic 
CRs other than hydrogen: helium, carbon and iron. The distribu- 
tion function (l44l i for the heavier elements is obtained by simply 
scaling the proton distribution function g(r, E) with the abun- 
dance with respect to the hydrogen;^fi of the corresponding ion / 
= H, He, C, Fe 



gi(r,E) ^Xig(r,E) 



(45) 



where E is the total kinetic energy of the CR ion. For the CR 
abundances we assume the standard values of the galact ic CRs: 
iLi.ArHe : A'c : ;rFe = 1 : 0.1 : 10"- : 10"^ (Maurin & TMiei 



i 



20031) . The low-energy cut-off for the heavier ions has been 
derived in the same way as for Galactic CRs, i.e. scaling the 
proton value ZiJ^j^ = 50 MeV by the ion atomic mass, E^^^^ 
=_My_£E^. but different parametrizations are also possible 
dNath & Biermann 19941) . 

The CR intensity I{E) (to which we usually refer as "spec- 
trum") is related to the distribution function by the following 
expression 



//(£) 



47r 



gi(rc, E) 



(46) 



where v - ficis, the velocity of the ion. The distribution function 
gi is calculated for a given distance rc from the center of the 
cluster. We adopt the value corresponding to the cooling radius 
of the cluster, which defines the region where the majority of the 
hadronically-generated y-ray luminosity originate. 

When the CR intensity is known, we can estimate the 
total lifetime against CR bombardment, tq. As an example, 
we calculate tq in two ne arby cooling flow clusters from the 
iPfrommer & EnfilinI (l2004h sample, A85 and Virgo, whose main 
properties are listed in Table |3] A85 is one of the farther cluster 
in the sample, with low electron density rif., large core radius rc 
and high central temperature Tq. Virgo is the closest cluster to 
us, with quite high electron density, very small core radius and 
hig h central temperature . 

iPfrommer & EnBlinl(l2004h calculated «cr,o for different val- 
ues of the proton spectral index a^. We adopt the value corre- 
sponding to ffp = ttinj + Q-ditf = 2.7, with Q'diif ~ 1/3, because is 
consistent with the spectrum of the CR sources (AGN + SNRs) 
and allows a direct comparison with Galactic CRs. 

The results of our calculation are shown in FiglT2] In the 
galactic clusters under examination the cosmic-ray density is en- 
hanced with respect to the solar neighborhood, resulting in life- 
times which are shorter by about 2-3 order of magnitude (cf. Fig. 
[8]). The longer time constants obtained for Virgo are the conse- 
quence of its lower CR density with respect to A85. 

In our previous work (MJTb) we showed that collisions with 
thermal electrons and protons represent an efficient destruction 
mechanism for PAHs immersed in an ionized gas of ~ 10'' K. 
Since the ICM is a moderately tenuous hot plasma, it would be 
not surprising to find a quite short lifetime thg for the PAHs 



Table 3. Relevant parameters and CR normalization factor «cr,o 
for the two cooling flow clusters A85 and Virgo from the 
Pfrommer & EnBlin (2004) sample. 



Cluster 


z 


{h-^ kpc) 




T 

(10' K) 


"CR,0 

(/^vo' cm-3) 


A85* 
Virgo^ 


0.0551 
0.0036 


45 
1.6 


3.08x10-^ 
1.5x10-' 


6.4 
1.4 


9.9x10"^ 
4.2x10-' 



(a): From lPfrommer & EnBlii] ||2()04|). w ith q-„ = 2.7. 

{b): Parameters from'Mohr et al.' ('1999") and'Oegerle & H illl j2()()ll) 
(c): Parameters from .Matsushita et al. ( 2002) and Ebeling etalJdl998h . 
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Fig. 12. Total PAH lifetime against CR bombardment in the cool- 
ing flow galaxy clusters A85 and Virgo, compared with the sur- 
vival time, Thg, in a hot gas. The value adopted for the lowest 
energy of the CR spectrum, Emin, is 50 MeV/nucleon (see text), 
otherwise as Fig. [8] 



dispersed in the medium, due to such collisions. We calculate the 
PAH lifetime in the hot intra-cluster medium of A85 and Virgo 
using Eq. 27 in MJTb and the temperatures and densities listed 
in Table [3] In both clusters, PAHs embedded in the hot gas are 
rapidly destroyed on timescales of 1 - 10^ yr In cold embedded 
gas, the lifetime of PAHs is set by CRs to be -10*" - 3x10** yr 
(Virgo) and -lO*' - 3x10' yr (A85). 

The existence of PAHs in the intra-cluster medium is not yet 
well constrained observationally. In fact to the best of our knowl- 
edge a detection has not been reported, but hopefully future ob- 
servations will be able to provide more insights. Many of the 
quantities involved in our calculation are not well constrained, 
in particular the exact amplitude, shape and composition of the 
CRs spectrum in the cluster and the value of the parameters £0 
and T(). Because the PAH lifetime in galaxy clusters set by CRs 
is still short, a replenishing mechanism would be required to ex- 
plain their observation in such environments. We face here an- 
other example of the well-known conundrum about the discrep- 
ancy between the PAH/dust lifetime and the injection time scale 
into the interstellar medium (2.5x10^ yr). 



8. Conclusions 

We have extensively investigated the stability of PAHs against 
CR ions (H, He, CNO and Fe-Co-Ni) and electron bombard- 
ment in both galactic and extragalactic environments. We con- 
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sider CR particles with energy between 5 MeV/nucleon and 10 
GeV. Collisions can lead to carbon atom ejection, with a conse- 
quent disruption and destruction of the molecule. The effects of 
CRs were then compared with the other destruction mechanisms 
discussed in our previous works: PAH processing by shocks with 
velocities between 50 and 200 km s"' (MJTa) and collisions with 
thermal ions and electrons in a hot gas (T = lO^* - 10'^ K, MJTb). 

An ionic collision consists of two simultaneous processes 
which can be treated separately: a binary collision between the 
projectile ion and a single atom in the target (nuclear interaction) 
and energy loss to the electron cloud of the molecule (electronic 
interaction). In the high-energy regime considered here the nu- 
clear stopping is totally negligible, and the energy loss process 
is dominated by the electronic interaction, well described by 
the Bethe-Bloch equation. The interaction of PAHs with high- 
energy electrons can be treated in term of a binary collision be- 
tween the incident electron and a single nucleus in the target. 

The CR spectra we adopt in the solar neighborhood are based 
on measurement near the Earth but corrected for the influence of 
the Heliosphere (solar modulation). To estimate the CR varia- 
tion across the disk and in the galactic halo we adopt specific 
models based on y-ray measurements. In external galaxies we 
scale the overall CR density with the star formation rate of the 
galaxy, adapting scale lengths and scale heights appropriate for 
the Milky Way. 

We find that the timescale for PAH destruction by CR ions 
depends on the electronic excitation energy Eq and the amount 
of energy available for dissociation. Small PAHs are destroyed 
faster, with He and the CNO group being the more effective pro- 
jectiles. CRs are able to process PAHs in diffuse clouds, where 
the destruction due to interstellar shocks is less efficient. For 
electron collisions, the lifetime is independent of the PAH size 
and varies with the threshold energy Tq. The minimum life- 
time is 1.2x10'^ yr, longer than the Hubble time. Such a long 
timescale excludes CR electrons as an important agent for PAH 
destruction. 

PAHs have been detected both in the halo of normal galax- 
ies like NGC 891 and in the outflows of starburst galaxies like 
M82. Our work shows that in both these environments the life- 
time against CR bombarding of large PAHs (A^c > 150-200) is 
comparable to or longer than the circulation timescale between 
disk and halo and the starburst lifetime (~200 Myr and ~ 20 Myr 
respectively). PAHs dispersed in the hot gas filling the galactic 
halo and the starburst outflow are rapidly destroyed by collisions 
with thermal ions and electrons, but this mechanism is inefficient 
if the molecules are isolated from this gas in denser cloudlets. 
CRs can access the denser clouds and together with X-rays will 
set the lifetime of those protected PAHs, which can be used as a 
'dye' for tracing the presence of cold entrained material. 

In cooling flow galaxy clusters like A85 and Virgo the cos- 
mic ray intensity is remarkably enhanced with respect to the so- 
lar neighborhoods, as a consequence the PAH lifetime is much 
shorter Nevertheless, the survival time against CR bombardment 
is at least two orders of magnitude longer than the PAH lifetime 
in a hot gas (1 - 10^ yr), implying that the molecules will be 
rapidly destroyed in the gas phase of the ICM. They could sur- 
vive if protected in some cold entrained material and in this case 
the PAH lifetime will be set by CRs. Future observations would 
hopefully provide more insights about the validity of our predic- 
tions. 

The major source of uncertainty in the determination of the 
time scale for PAH processing by CRs resides in the choice of 
the nuclear threshold energy To (for electron collisions) and the 
fragment binding energy £0 (for ion collisions). Our conclusions 



are robust despite the large variability in the PAH lifetime in- 
duced by the incertitude on the above parameters. Nevertheless 
this variation emphasizes again the importance of a better deter- 
mination of these quantities. 

We find that thermal ions/electrons in a hot gas are much 
more effective in destroying PAHs than CRs. This is due to the 
fact that the stopping power of the thermal ions/electrons un- 
der consideration {T ~ 10^ K, «e ~ 0.1 cm"-') is high, close to 
its maximum value, and allows the transfer into the molecule 
of enough energy to have the dissociation probability close to 
one (the dissociation probability increases with the transferred 
energy). Almost any ion/electron is able to destroy a PAH, and 
the resulting destruction efficiency is very high (cf. MJTa and 
MJTb). 

On the other hand, in the energy range we consider here 
for CR ions (5 MeV/nucleon - 10 GeV), the stopping power, 
and then the energy transferred into the PAH and the dissoci- 
ation probability for PAHs, decrease rapidly for increasing en- 
ergy of the incoming ion. For CR electrons (5 MeV - 10 GeV) 
the cross section for carbon atom removal is almost constant but 
very small (~ 6x10"^'' cm^). For both ions and electrons the CR 
spectra are decreasing functions of the energy. The combination 
of these two factors implies that only few CRs are able to destroy 
PAHs, resulting in a destruction efficiency very low compared to 
thermal ions/electrons in a hot gas. 

Our results show that CRs set a timescale for the destruction 
of PAHs with less than a ^ 100 C-atoms of only 100 Miflion 
years in the interstellar medium. Larger PAHs or very, very small 
dust grains - with A^c > 100 atoms - are predominantly pro- 
cessed by interstellar shock waves (rather than CRs) on a very 
similar timescale (Tshock - 150 Myr; MJTa). 

PAH molecules are modelled to be efficiently formed as 
molecular intermediaries or sideproducts of the soot formation 
process in the st ellar ejecta from C-rich objects, in particular C- 
rich AGB st ars dFj enklach & Feigelson 1 9891: ICherchneff et aP 
1 19921 120001: ICherchneff .2000.: .Cau 2002i) . If we assume fliat 
AGB stars are the primary source for PAH replenish ment in the 
ISM, the injection timescale, rfoimation - 3x10^ Myr dJones et al.l 
1994) is much longer than the lifetime estimated above. Thus, 
there is a factor of about 50 discrepancy between the injection 
and destruction timescales for PAHs in the ISM. This parallels 
the discrepancy betwe en the iniection an d destruction timescales 
for interstellar grains (I Jones et al. As for dust grains, the 

conclusion seems to be inescapable that PAHs must be reformed 
rather efficiently in t he ISM itself . 

As suggested by iJones et al ] (ll996h . interstellar PAHs may 
be the fragmentation products of carbonaceous grain collisions 
in shocks. Experiments mimicking the effects of grain-grain col- 
lisions support the formation of PAHs and fuUerenes as well 
as carbon chains and clusters from Hyd rogenated Amorphous 
Carbon (HAC) grains (IScott et al.lll997h . This would replenish 
the PAH population on a similar timescale (^ 100 Myr) as the 
destruction timescale. With the numbers quoted above, if all of 
the carbon were injected into the ISM in the form of HAC grains, 
this scheme would (almost) be in steady state with 5% of the el- 
emental C in the form of PAHs. However, this comes at the ex- 
pense of the rapid destru ction of carbon grains, as s eems to be 
suggested by the work of ISerra Dfaz-Cano & JonesI (l2008h . So, 
this would have to be fit into a general model where gaseous car- 
bon accretes and reacts with carbonaceous (e.g., HAC ) grains in 
intersteflar clouds (cf., lJoneslll990l:lJones et al.lll990l) . 

This growth process is then balanced by grain-grain colli- 
sions in strong shock waves in the intercloud medium produc- 
ing PAHs that are then destroyed through CR interaction in the 
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diffuse ISM. While such a model could work numerically, the 
basic premise — rapid accretion and reaction of gaseous car- 
bon on grains forming predominantly HAC-like layers - has not 
been demonstrated in the laboratory. Moreover, observationally, 
depletion studies suggest that carbon does not seem to partake 
in the rapid ex change between gas and solid phases in the ISM 
(lTielensll2009l) . Alternatively, rather than interstellar PAHs re- 
sulting from chemical processes in stellar ejecta - either directly 
or through the shattering of carbon grains - these species may be 
largely formed from smaller carbon molecules through chemical 
processes in the ISM. 

Theoretical studies show that ion-molecule reactions can 
form benzene-like species in the dense and warm environments 
of pre-Planetary Nebulae (such a s CRL 618) and in proto- 
planetary disks around young stars ( Herpin & Cernicharoll2000l : 



IWoods et al.ll2005tlWoods & Wmacvl2007ir However, such spe- 
cial environments are few and far apart and not important in a 
model for the global presence of PAHs in the ISM. Models pre- 
dicts that ion-molecule chemistry in cold dense molecular clouds 
can produce benzene but at peak abundances at early times of 
10"^. Likely more relevant are the stead y state abundances a t 
long times and those are much less (10'"^: iMcEwan et al.ll999h . 
So, the gas phase formation of PAHs in molecular clouds does 
not seem a very promising route. Essentially, the basic building 
block of PAHs (acetylene) is not very abundant and there isn't 
much free carbon around to drive the chemistry since most of it 
is locked up in CO in molecular clouds. In conclusion, the ori- 
gin of interstellar PAHs is still clouded in quite some intriguing 
mystery. 
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